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An approximation to General Relativity is presented that agrees with the Einstein field equations 
up to and including the fourth post-Newtonian (PN) order. This approximation is formulated in 
a fully constrained scheme: all involved equations are explicitly elliptic except the wave equation 
that describes the two independent degrees of freedom of the gravitational field. The formalism 
covers naturally the conformal-flat-condition (CFC) approach by Isenberg, Wilson, and Mathews 
and the improved second PN-order exact approach CFC+. For stationary configurations, like Kerr 
black holes, agreement with General Relativity is achieved even through 5PN order. In addition, 
a particularly interesting 2PN-exact waveless approximation is analyzed in detail, which results 
from imposing more restrictive conditions. The proposed scheme can be considered as a further 
development on the waveless approach suggested by Schafer and Gopakumar [Phys. Rev. D 69, 
Q ■ 021501 (2004)]. 
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I. INTRODUCTION 

' The solution of the Einstein field equations for inspiraling compact binaries is one of the great challenges of Numerical 
Q-< Relativity, e.g., [If. To fully succeed with this problem realistic initial data are needed which do not include incoming 
(spurious, "junk") gravitational waves. If one could follow up the evolution from the very beginning of the inspiraling 
process, i.e., right from the infinite separation of the both objects and without gravitational waves present at this 
initial stage, the configuration to start with at any finite instant of time would be well under control just by calculating 
the history up to that instant of time. Lack of this piece of information requires the construction of initial value data 
^ . by other means. 

' On the other hand, there are several astrophysically relevant scenarios, like supernovae explosions or initial stages of 
, a binary black hole, which can be described very accurately with approximated versions of Einstein equations. These 
■ two issues (initial data generation and construction of truncated schemes of General Relativity) are deeply connected 
since by extending to all times the conditions imposed on the freely specifiable part of the initial data set, one can 
£N) ■ obtain an approximated version of Einstein equations. 

As is well known, in the initial slice, we are free to choose four gravitational degrees of freedom from the initial data 
set (the two independent field degrees of freedom and their time derivatives or, their canonically conjugate momenta), 
whereas the remainder gravitational objects must be obtained through the resolution of the constraint equations and 
J> , coordinate conditions. Resorting to a canonical formulation of General Relativity, e.g., the Arnowitt-Deser-Misner 
(ADM) one @ - to be used throughout in this paper -, these independent gravitational degrees of freedom get clearly 
rS \ identified. The four degrees of freedom are encoded in the two independent metric-field components hJ T and their 

?H 3 

canonical conjugate 1 g C 7r g 7r TT' where TT means transverse and traceless with respect to an auxiliary 3-dimensional 
flat metric 5ij with i,j — 1, 2, 3. Therefore, the initial value problem reduces to the fixation of the spatial dependence 
of these two objects at a finite initial time. 

In the Isenberg- Wilson-Mathews conformal-flat condition (CFC) approach to General Relativity 0,0], hj^ and 
its time derivative hj^ are chosen to be vanishing. Hence, this scheme is very appropriate to construct conformally 
flat initial data. On the other side, understood as an approximation to General Relativity so that the mentioned 
conditions (hj^ — hj^ — 0) are valid for all times, this formalism reproduces the evolution of the gravitating system 
at 1PN order, which implies deviations from General Relativity at 2PN. 

Even though, near-zone post-Newtonian (PN) calculations show that for inspiraling binaries with no-incoming 
radiation, hj^ can never be exactly zero, see, e.g., but contain conservative and dissipative terms which destroy 

CFC. The first piece of deviation from CFC appears at 2PN order and is of conservative type. The higher-order pieces 
appearing in the cited papers are 2.5PN (dissipative), 3PN (conservative), and 3.5PN (dissipative). Therefore, many 
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attempts have been made to construct conformally nonflat schemes, see for instance [10|, In particular, regarding 
works done based on PN corrections, the so-called CFC+ approach, as detailed in [12], generalizes the CFC framework 
by just correcting for the additional 2PN piece. In the scheme [l3[ also 2.5PN terms have been considered and were 
recently further developed in [H], The 2.5PN and 3.5PN pieces have been added within the CFC-type skeleton 
approach [r| in [l7j . However, in this last reference the radiation reaction has been taken into account even through 
6PN order. Herein, the 4PN, 5PN, and 6PN orders are both conservative and dissipative, where the dissipative parts 
result from the peculiar tail structure of the radiation process. 

The separation in dissipative and conservative terms allows the separate calculation of the dissipative and conser- 
vative parts of the metric coefficients. The aim of the waveless approximations to General Relativity beyond CFC is 
the determination of conservative terms in the gravitational radiation process. Related with waveless approaches are 
the near-zone helically symmetric schemes introduced by [l8j , which also do not allow for emission of waves and have 
been intensively studied in the recent past, e^., [l9l - l2l| . The fully constraint formulation [22J has found deeper inves- 
tigations and waveless applications in (23l . |24| . A comparison between helically symmetric and waveless description 
of binary systems has been performed in, e.g., [25j . 

In the present paper, a modification of the waveless approach by [26] is developed which supplies a treatment of 
radiating systems in agreement with full General Relativity through 4PN order and even up to 5PN order under 
nonradiating conditions. The CFC and CFC+ approaches are well covered as well as treatments including outgoing 
radiation [l3l - [T5l . [27j. This approximation is formulated in a fully constrained scheme. All involved equations are 
explicitly elliptic except the wave equation obeyed by hj^ that describes the two degrees of freedom of the gravitational 
field. This hyperbolic equation must be fed with freely specifiable initial data hJj T and hJ^Q. 

In addition, we analyze in detail a particularly interesting 2PN-exact waveless approximation to General Relativity, 
which is a particular case of our 4PN approach. More precisely, it results from imposing the "conjugate" conditions 
to CFC, namely 7r^ T = 7r^? T = 0. 

Under stationary conditions only integer PN orders are different from zero. We will keep to this integer PN-counting 
also under nonwavcless conditions where half-integer orders appear. This particularly means that our proposed 4PN- 
exact approach is, in precise terms, a 4.5PN-exact approach as well. 

The article is organized as follows. In Sec. |TT] the Einstein field equations are shown in the context of the ADM 
formalism. Section IIIII presents explicit elliptic equations that are valid in full General Relativity and allows one to 
solve for all the geometric objects but the transverse and traceless part of the metric and of its conjugate momentum. 
In Sec. IIVI the transverse and traceless part of the conjugate momentum is assumed to be vanishing, which leads to 
a 2PN-exact waveless approximation scheme. Section [V] generalizes the previous scheme to a 4PN-exact framework 
suggesting another 2PN-exact waveless approach. We show that this formalism reduces to the well-known CFC+ 
approach under the corresponding conditions in Sec. IVII The particular matter model of a perfect fluid is presented 
in Sec. IVIII Finally the conclusions are drawn in Sec. IVIIII 

II. EINSTEIN FIELD EQUATIONS IN ADM FORMALISM 

In the ADM formalism of General Relativity, the spacetime line element is split into the (3 + 1) form, 

ds 2 = -a 2 c 2 dt 2 + jijidx 1 + /3 i cdt)(dx j + P J cdt), (1) 

where a is the lapse function, fi l the shift vector, 7^ the induced metric on a three-dimensional spatial slice 

parameterized by the time coordinate t, and c is the speed of light. The three-metric 7^ and its canonical conjugate 
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16 c ^ G ir v , which is a contravariant symmetric tensor density of weight +1, satisfy the Hamiltonian and momentum 
constraints @], 

^ 1/2R = ^vf^^ - **) + -7 1/2 « 2 r 00 , (2) 

(-27^ =) - 2tt^ + ir kl lk u = K^^aT , , (3) 

where R is the curvature scalar of 7 the determinant of 7^, 7T*j = 7jfc7r lfe , and tt = 7rV On the other hand, 

T 00 and T°i are the components of the unspecified four-dimensional stress-energy tensor for the matter source T^ v . 
The canonical conjugate 7r y is related to Kij, the extrinsic curvature of by tt^ = —^ 1 / 2 {^ tl ^i m — 'y l i^ lm )Ki m , 
where 7 y is the inverse metric of 7^. In the above equations, a partial derivative is denoted by a comma, whereas 
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stands for the three-dimensional covariant derivative, and k = \a , with G the Newtonian gravitational constant. 
The three- metric and its canonical conjugate evolve in accordance with the following evolution Eqs. [2|], 

tt 1 \ = -ia 7 1/2 (2R lJ - 7 lJ R) + ia 7 ~ 1/2 7 y '(27r M 7r fei - tt 2 ) - aj- 1/2 {2TT lk ^ k - irir 11 ) 

+ 7 V2 (a l« _ yi a|fc fe) + _ _ n ik^ + I a ^/^ 7 ji Tku (4) 

and 

lijfi = a~f- 1/2 (2irij - 7ij7r) + + (5) 

where Ry is the Ricci tensor associated with In this paper, we raise and lower indices of three-dimensional 

objects with 7 IJ and jij respectively. 

The ADM coordinate conditions, which generalize the isotropic Schwarzschild metric, read 

tt" = 0, (6) 

where, and from here onwards, repeated covariant or contravariant indices imply the usage of Einstein summation 
convention, and 

Hi = *&ij + h??, (7) 

where \& is a conformal scalar and hj? the transverse and traceless part of the three-metric jij with respect to 
the Euclidean 3-metric By definition, hf^ satisfies hj^ = hj^ = 0. The conformally-flat condition hj^ = 

gives a simple expression for the three-dimensional curvature scalar, 7 1 / 2 R = —8 A^ 1 / 4 , where A stands for the 
three-dimensional Euclidean Laplacian. The gauge fixing equation for 7^ ([7]) can be rewritten in a differential way, 

37«,3 ~ 733,; = 0. (8) 

Equation © results in the covariant trace of 7r y of the form tt — ir^hj^ . Taking into account the space- asymptotic 
properties tt % ^ ~ I/r 2 and hj^ ~ 1/r, the gauge condition ^ turns out to mean asymptotic maximal slicing 
(K - 1/r 3 ). The gauge conditions © and 0, or ©, are very close to the well-known Dirac gauge conditions which 
include maximal slicing. In fact, condition is identical with the corresponding Dirac gauge condition to linear 
order in 7^ — 5ij . 

As will be made explicit in the next section, the function ty, with — 1 ~ 1/r at asymptotic infinity, and the 
longitudinal part of the momentum 7r y are determined using the Hamiltonian and momentum constraints, given in 
Eqs. © and ©, by elliptic equations. On the other hand, the elliptic equations for the lapse a, with a — 1 ~ 1/r, 
and the shift f3 l , with f3 % ~ 1/r, result from the evolution Eqs. ((4]) and ([5]), respectively, after applying the coordinate 
conditions ([5]) and Since we want to provide a set of equations that can be numerically solved, we also will take 
care that all the sources decay at least as 1/r 4 at asymptotic infinity, avoiding in this way convergence problems. 

Finally, it is important to stress that our variables are those dictated by a canonical formulation. This fact clearly 
makes them form a fully consistent set of independent variables. Of course, other sets of variables may also be 
chosen, like the one presented in Ref. [22j or the Kol-Smolkin one [28j |. which is connected with the Landau-Lifshitz 
decomposition of the metric [2!| and has extensively been used in, e.g., [30| . Even though, it should be pointed out 
that only within a radiation-type ("Coulomb- type") set of coordinate conditions those variables fulfill a mixed set of 
elliptic and evolutionary (hyperbolic) field equations instead, e.g., a purely hyperbolic system. 

III. EXPLICIT ELLIPTIC EQUATIONS 

In this section we will present explicit elliptic equations to solve for all the objects but the transverse and traceless 
parts of the three-metric 7y and of its corresponding momentum 7r u . No assumptions will be made, hence these 
equations will be valid in full General Relativity. 

A. Equation for the shift 

Taking the time derivative of condition ([5} and making use of Eq. (O, we obtain an elliptic equation for the shift, 

(/%),3 + (%),3 " §(%),< = -(a^ 1/2 (2n lJ - 7y 7r))j + ^(a7~ 1/2 (2^ - 7^)),^ (9) 
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It can be more explicitly written down in the following way, 

Aft + \Pj.ji = -l(l jk ^T kll ),i + 2{ r f k p i T m ),x + (crT 1/2 7^)d + ^("7- 1/2 (2^ - ljj7 r)), 



2^(«7- 1/2 ),i-^ 



(10) 



where we have applied the momentum constraint §3§ to replace the term %%j j in order to guarantee an adequate 
fall-off behavior (~ l/?" 4 ) at asymptotic infinity of the right-hand side of this equation. Here one should use the usual 
definition for Christoffel symbols, 



and the exact relations 



77' 



Fijk = ^ (jij.k + Jki.j - Jjk.i), 



\'\, +Hij, 



(11) 

(12) 
(13) 



being x = 'J 2 — ^hf^hj^ and = hJ k T hJj^ — ^hj^, for the inverse and the determinant of the metric, respectively. 



B. Equation for the lapse 

Combining relations (U) and © it is straightforward to obtain an elliptic equation for a, 



rv i J 



7 a 



= -aH +(7 dik — 2'y J cry ir J ir jk - —7 7r - 27 ' 7T p , 



(14) 



where we have made use of the Hamiltonian constraint © to get rid of the Ricci scalar. Note that the last term of 
the first line is corrected from Ref. (26|. The left-hand side of this equation can be given as, 

(y™ 7 « _ yy ) ( adm - 7 fci r ym a, fc ) , (15) 

where, the inverse metric combination can be rewritten making use of the formulas presented in the previous section, 

1 



2y* 2 

(7 7 -7 7 ) = 2~ <>jm t —~ 



Therefore, the equation for the lapse takes the explicit form 

1 



(x + Ha) Hj m — , 



(16) 



Aa 



2 



(X + Ha) Hj m — HijHi, 



k u 



7 



2x*" 

K 



-aR + (7 Oife - 27 J a 7 7T J Wj k - —7 7T - 27 7 7T ^ 



+ £ a ( 7 V n Tzm + 7"a 2 T °) 



(17) 



The only term that is left to be known in this equation in terms of metric components is the noncovariant trace of 
the Ricci tensor R M . This tensor can be written in the following way, 



2 7 R 4 



H k l(j k l,ij + Jij.kl — Jkj.il — Jil.kj) 



+2^^ kl Y LP (^nilFp k j ~ FnijTpkl) , 



(18) 
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and its trace is then given by, 
R" = 



2x 3 



A* 



X 

2 7 3 



H i7n^l jrn 



(AhJ? + SijA* + * 



'2l 



J tm J 3m HklfaklM + lij,kl — Jkj,il - lil,kj) 



+ i im i jm i kl i np ^mir pkl -r mj r pkl ). 



(19) 



In this expression one should replace the Laplacian of the conformal factor that appears in the first term of the 
right-hand side with its corresponding Eq. pip, which will be presented in the next subsection, in order to enforce a 
decay rate of ~ 1/r of the right-hand side of Eq. (fl7|). 



C. Equation for the conformal factor 



We will use the Hamiltonian constraint © to solve for the conformal factor ^. The Ricci scalar can be written in 
the following way, 

R = (/y <*yj _ 7 V)(tf )tf a w + hjfia) + i kl i np ^{T ml T pk] - T nij T pkl ), (20) 

where the first- and second-order derivatives of metric components are clearly separated. Expanding it we get the 
following more explicit expression, 



(x + H kk )Hij ~ H ik H jk 



1 



Hi k H,ji — HuH, 



ij-n k i 



hj? u + i kl i np f\r ml r pkJ - r nij r pkl ) 



Now it is straightforward to introduce this expression in the Hamiltonian constraint obtaining, in this way, an elliptic 
equation to solve for the conformal factor, 



A* = KtHiiAhJ? 

2* 2 3 13 2x* 2 



(x + H kk )Hij — H ik Hji 



1 



2x* 2 

7 

4x* 2 



Hi k Hji — HijH, 
(2^%^- - ir 2 ) - 



1,-1 



n ij, k l 



2 X * 



r7 H 7 n V(r n «r pfej -T nij T pkl ) 



r a 2 7 2 T u 



(21) 



D. Equation for the longitudinal field momentum 



Because of the gauge condition ©, we can perform the following decomposition for 7r u , 



'TT' 



where the longitudinal part tt ij can be given in terms of a vector field, 



2 a r 



An elliptic equation for this vector can be obtained from Eq. ([3]), since n lJ t j = An' 

1 , 



3 " i*J ' 



Att 1 



'1 1 kjlir - 2 7 «7 J j 



(22) 



(23) 



(24) 



In summary, we have fixed the eight (out of 16) degrees of freedom in the initial slice through gauge conditions 
([6|) and {2} together with the constraint equations, which are implemented by the elliptic equations for the conformal 
factor (f2"Tj) and for the longitudinal part of the momentum (IM1) . Finally, Eqs. (fTU|) and (TIT)) for the shift and lapse 
guarantee that the mentioned gauge conditions are conserved through the evolution. 

The only objects that are then left to be fixed are the transverse and traceless part of the metric hJ 3 T and its 

conjugate momentum 7i%j? T . The four components of these tensors describe, in this setting, the two physical degrees 
of freedom of the theory. This means that they are freely specifiable on the initial slice but the conditions imposed 
on them will restrict the physical problem we are dealing with. 
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IV. A 2PN-EXACT WAVELESS CONDITION: yriL = 



Thus, in order to close the set of equations and obtain an approximated version of Einstein theory, we will introduce 
the following restrictions, 

7r^! T = and thus, 7r^ T! o = 0. (25) 

This assumption is exact at 2PN order since, in the general scenario, 7nwr is vanishing up to 0(c ) and it can be 
understood as the conjugate condition to that of the CFC case, where /i^ T and its time derivative are assumed to 
be vanishing. Therefore, in this approach, the conjugate momentum ir lJ will be purely longitudinal and, hence, the 
vector field 7t z will contain all the needed information in order to reconstruct it (|23[) . Thereby now the only missing 
equation is that for hj? and it will be obtained from the TT part of Eq. Let us denote by A ij the right-hand 
side of that equation, 



A ij ee -aj 1/2 R lJ + Y ij 



where we have defined, 



Y l i EE _^(2 7 «7r W 7r fe i - Air ik ir' k + 2^ - n 2 j ij ) 

+ 7 (7 7 ~7 J 7 )(afem-7 Infcro",;)-",; 



+ (ir ij /3 m ) tTn - Tr mj P\ m - 7r mi p\ m + - ai 1/2 {f l ^ m T lm + ^a 2 T 00 ), (26) 

and use has been made of the Hamiltonian constraint ([2]) in order to remove the Ricci scalar from this expression. 
The term a,^ has been added and subtracted from the previous definition in order to ensure a decay rate of 1/r 4 for 
Y l i at spacelike infinity. 

On the other hand, raising indices in formula (|18[) . the contravariant Ricci tensor is given by, 

R« = (Ah?? + 9, a + fyA*) + XV, (27) 



with 



X ij ee —^L [AhJ^ + 4mA* + * )fcm ] Ix^miHjk + xSkjH im + H im H jk 



2^7 7^ (ylnpikin ~t~ 7km, np ^nm^kp 7fcp,nm)-^np 
_|_ J lk J : ' m Y S ^ nP {T n ks^prm - r„fe m r prs ). (28) 

With these definitions at hand, one can write down as, 

1 f a X 3 ,\ /a ,tt 



1 



-(A/^ + ^+^A*), (29) 

where, again, terms in the second line have been added and subtracted for future convenience. Because of the gauge 
condition ([5]), A 13 is tracefree on shell. Even though, in order to make an eventual numerical evolution of this system 
more stable, we will make it explicitly tracefree. After that we will need to obtain the transverse part of the resulting 
object. In order to do so, we take into account that full-derivative terms do not contribute to the transverse and 
traceless part, thus both a ^ and * can be eliminated, whereas explicit pure-trace terms (those proportional to <5y) 
can be included in the trace part A. Therefore, the evolution equation for 7r^ T ((4]) is now rewritten as, 

4 T ,o = i(A^T_ jB ^ ) (30) 

where the source term is given by, 

B% T ee B ij - l^B - V\i - V\j + , (31) 
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being B the noncovariant trace of B^ , i.e., B = B n , and 




Finally, the vector field V' L encodes the longitudinal part of B % i and hence obeys the following equation, 

W + l Vj ,ji = Bij ,i ~ \ B ,i ■ ( 33 ) 

Imposing now our assumption (TrSL = tt^ t q = 0) on Eq. (|3H)) . it is straightforward to get an elliptic equation for 
'hj > 

Ah$?=B i 4 T . (34) 

This assumption must also apply on other elliptic equations presented in last section, obtaining in this way a closed 
system of equations which does not contain 7ivj? T . Because of the terms that have been added and subtracted previously, 
the source B^ T decays properly (~ 1/f 4 ) for an eventual numerical resolution of this equation. In particular, note 
that the terms contained inside the first brackets of the definition of _B iJ (|32[) have been constructed such that their 
combination decays as ~ 1/r at asymptotic infinity and, moreover, it is of order 0(c~ 2 ), as can be verified using the 
expressions (|3"5|) below. 

In summary, the gravitational metric functions are encoded in the objects {j3i, a, hj^} which are connected with 

the auxiliary variables {7T 1 , V 1 } and the corresponding elliptic equations are, once imposed the condition ir^ T = in 
all of them, ([ID ]) . ([T7 |) . ()53 ]> . and (|55 ]) , respectively. 

This set of equations has two main applications. On the one hand, it can be used to obtain initial data, such that 
= 0, for a subsequent numerical evolution of the full Einstein equations. On the other hand, it could be 
solved as a self-contained theory of gravitation, which would not have waves, but would approximately describe the 
motion of the astrophysical bodies. 

Regarding this last application, we would like to end this section by analyzing the accuracy of our theory, that 
is, which is the error one would commit when solving our approximated set of equations with respect to solving the 
full Einstein equations. Since there is no approximation in the equations for the rest of the objects, the only error 
is introduced in the system via the assumption on n^ T . Looking at their respective equations, and recalling their 
leading order post-Newtonian series, namely, 

a = 1 + O (c- 2 ) , Pi = O (c~ 3 ) , tt« = O (c- 3 ) , 

* = l + 0(c- 2 ), l^ = 0(c-i), (35) 

one can check that if 7r^ T ( n ) were computed properly, where the subscript (n) denotes the coefficient of 1/c™ in 
the corresponding post-Newtonian expansion, the metric components {c»!( n +3), /3(„ + 2)> 7«i(n+i)} would be obtained 
correctly. In fact the conformal factor is calculated up to ^( n +3)i otherwise we could not compute the lapse up to the 
mentioned order. As we have already commented, the assumption we have proposed is correct up to 2PN order, i.e., 
we are computing 7r™ correctly thus, by using the proposed scheme, {c*(6) j /?(5)> 7y (4)} can be properly determined. 



V. A 4PN-EXACT APPROXIMATION 



In this section we will generalize the scheme we have presented in the previous section, which will give rise to a 
formalism that agrees with General Relativity up to 4PN order. Contrary to our 2PN scheme, this more general 
approach will indeed contain gravitational waves since the main equation to obtain hj^ will be of hyperbolic nature. 

Let us write the evolution equation for the spatial metric ([5]) in the form 

Tij',0 = Cij , (36) 

where 

dj = a 7 - 1 / 2 (2 7 r lJ - 7 ^7t) + /% + %. (37) 

Using the gauge conditions on 7^ , that are given by Eq. |(7) , we obtain the evolution equation for the transverse and 
traceless part of the metric, 

hjjfi = C ij - 7 (38) 
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where C u = 3r](C il - Y j J tk h]^ ) with C u = -cry" 1 /^"*' + 20 and 77 = This equation can then be 

rearranged so that ir lJ follows in terms of hj^ in the form 



7 V2 

~2^ 



■ / im in _ ik jk Ira ZtixtTT 

+ (7 7^ -i)7 f 7 7 )h m 



(39) 



where the matrix M lJ is given by 



TT 



(40) 



Remarkably, the matrix M tJ m deviates from the unit identity matrix in the quadratic order of hJ T only [26| . Taking 
into account that hj? is of 2PN order (f3"5j). neglecting the quadratic terms of /i? T in M lJ \i is a truncation at 5PN 
order only. Thus, we can write down 



•yl/2 

~2o~ 



+ (7 7 J — 777 7 J 7 7 



ran.O 



(41) 



which is a 4.5PN-exact relation. Keeping only the leading order of the time derivative of hj? and calculating the 
transverse and tracefree part of this relation, we obtain, 



1 



*% r = ~h'g + D i t-DZ + O(c 



2 Hfi 



where 



and 



hold with 



D 13 = -- 



(42) 



(43) 



(44) 



(45) 



Note that the decay rate of the right-hand side of the last equation is 1/r 3 at asymptotic infinity. As we have 
commented, this could give problems when solving this equation numerically. Even though, since the term in question 
is a full divergence, the vector W l can be redefined including in it the terms that goes like ~ 1/r 3 , so that the source 
for the equation for the new vector decays properly. 

In order to check that Eq. (|42|) is valid up to 0(c~ 9 ) one has to take into account the fact that a 2n ^ n = l + 0(c~ 4 ), 
as can be verified with the first terms of their PN expansions: = l + 2U/c 2 and a = 1 — U/c 2 , U being the Newtonian 
potential. 



The (nonapproximated) evolution equation for ir^ T (|30p and the (approximated) one for hj^ (|42p . 



ft™ =2^-2^-^ 



(46) 



define our system of truncated evolution equations of the Einstein theory put into canonical form. Combining the 
mentioned two equations, and thus making an implicit change from a Hamiltonian to a Routhian framework (see, 
e -g-j @)) results into a second-order hyperbolic equation for /i^ T , 



(47) 
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which describes the propagation of the two gravitational degrees of freedom. Defining the source S 1 ? 



2 ( D ij - D l l J +B% T , the no-incoming radiation formal solution of this equation is given by the standard retarded 



with the retarded time i rot — t — \x — y\/c; for another representation of the retarded solution see, e.g., }31| . 

In summary, the system of equations to be solved is composed by the elliptic Eqs. (flO|) . ([17]) . ([2"T]) . ([24| . (f33|) and 
the wave Eq. (|47|) for ft-ij T . This set of equations can be understood as a truncated version of Einstein equations 
so, by assuming certain initial data for the hyperbolic equation, they provide an approximated solution of General 
Relativity. On the other hand, they can also be used to obtain initial data for a subsequent evolution via full Einstein 
equations. Here we propose an iterative scheme to construct such initial data for the metric coefficients. As a first 
step, one imposes the assumptions ir^ T = and hJ 3 T = in all elliptic equations and solves them including the 
equations for the matter variables for all past times. Then, with this information, one can solve the hyperbolic Eq. 
(H7J for hj^ via the retarded integral (and hence without any "instant-of-time" initial conditions) dropping the time 
derivatives of the right-hand side since they are of 0(c~ 6 ). In a second step, by solving again the elliptic system with 
the previously computed hj^ , and still 7r^ T = 0, a 2.5PN-exact solution {ci!(7),/9jL,7y(5)} is obtained. Note that this 
also 2PN-exact solution is essentially different from that obtained by solving the scheme proposed in Sec. IIVI This 
2.5PN-exact solution can now be used to compute k^tIS) properly through its truncated definition (|42l) . With this 
information at hand, we can calculate the sources (all the right-hand sides) of our system of equations to a PN-level 
of precision higher than in the second step. Therefore, by solving this system, a 3PN (and thus also 3.5PN)-exact 
solution can be obtained, {cy(g), 1% Repeating again this process, the metric components can be computed 
up to 4PN-exact order, {a(io), ffl g y 7ij(s)}- The limit of this procedure is given by the order up to which Eq. (|4"2"]l is 
valid. This iterative scheme generalizes a procedure which is under development in [lil [l5|. 

At this stage some remarks are in order. 

In this second-order formalism, ir^ T has to be replaced by hj? in all our equations using its truncated definition 
(|4"2"]l . This time derivative, as well as the first-order time derivatives on the right-hand side of the wave Eq. (|4"Tj) . have 
to be evaluated applying matter evolution equations, e.g., those that will be presented in Sec. IVIII 

Regarding the decay rates at asymptotic infinity, from Eq. (51) in Ref. [8|, it follows that at 3PN level the relation 
D lJ — Z?l = ( x I/(2)7r(3- ) ) TT is valid. Since for large radii fr^ includes the total linear momentum, this implies that the 

difference (D 1 ^ — ),o, that is present in (|47p. decays as ~ 1/r 4 at asymptotic infinity due to conservation of linear 
momentum. With respect to the implicitly present hJ } T in many sources, the conservative motion which enters into 
the source terms of the hyperbolic equation does not allow for a solution with time derivative decaying at least as 1/r 2 . 
Only at higher PN orders, when radiation damping is entering, the matter motion produces decaying wave amplitudes 
towards spacclikc infinity. In order to achieve a proper decay on our level of approximation a technical trick, similar 
to adiabatic damping, may be introduced as follows. All partial time and space derivatives c* M (with /i = 0, i) of hJ 3 T , 
that appear in source terms, should be replaced through exp(— er/ro)d fl , where r is the radial coordinate originating 
from the center of energy, ro a typical radial extension of the matter distribution, and e a small positive constant, 
< e << 1. At the end of the day, the final result should depend neither on e nor on ro. For a related treatment 
of retarded solutions see, e.g., [32j]. In conclusion, applying this procedure, through 3PN order, and even 3.5PN, all 
source terms in the elliptic and hyperbolic equations to be solved decay nicely at spacelike infinity. On the other 
hand, at the 4PN level tail terms will have to be controlled as it is evident from the Introduction. 

Finally, the Eq. (|41[) shows that for stationary configurations, like for Kerr black holes, where hj^ Q vanishes, our 
approach is 5PN exact, that is, it provides {a(ia), fi\ixy Jij (10)} accurately. This context suggests another 2PN-exact 
waveless approach through dropping hj^ ' in Eq. (|4ip and also all terms in Eq. (|47|) with time derivatives, which is 

equivalent to the assumption 7r^? T ,o = in Eq. (|30[) . In this way, a quasistationary approximation scheme is achieved 
which is 5PN exact for stationary configurations and covers both the lPN-exact CFC and the 2PN-exact CFC+ 
approach (see next section) . This latter approach is closest to (26( , where 7r^ T is also put equal to zero in Eq. (|30|) 
but Eq. (f39|) is kept exact. 




integral (see below for an iterative solution), 




(48) 



VI. RELATION TO 2PN-EXACT CFC+ APPROACH 



In this section we want to reduce our scheme to CFC+ order in order to compare with previous results in the 
literature. As has already been explained, the CFC case is nothing but allowing hJ T to be vanishing. This is the 
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only approximation that is done in this approach, so all the equations that we will present in this section will be fully 
correct in the conformal flat case just by removing terms with hj^ . On the other hand, the CFC+ approach considers 
also linear terms in hj^ hence, in this case, the following definitions apply, 



7 *i = - hf?, 7 = * 3 , X = * 2 , and n=\- (49) 

The time derivative of hJ T in Eq. (|41[) can also be neglected, including however a nonvanishing 7ivj? T which makes it 
a 2PN-exact approach different from Sec. IV, so for ir^ only 



7 1/2 / I- -i 



2a 
^1/2 

2a 



(50) 



is needed. In fact, this equation is only used in the CFC case since in the CFC+ scheme, it is only valid up to 1PN 
order c -3 counting in terms of integer PN orders. 
At this linearized order, our equation for the shift, 

Aft + = -^(7 J ' fc &r fe «),i + 2(y fc &r fciJ ), l - 2{a 1 ~ 1 l 2 Tr lJ ) j, (51) 

is trivially obeyed when taking into account last relation (|50|) . Hence, in order to obtain a meaningful equation for 
the shift, we introduce relation (|50p in the momentum constraint and obtain, 



A/? 4 + ~F d i = -2* (^) ^ tt" + Wt*TV (52) 
In the equations for the lapse (fTTJ) and the conformal factor (|21[) . the terms that survive at this order are given by, 



2 



Q 

7* ' 



(53) 



Aa = -^ T a„- +7 fc 'r ijJ a fe + g ^ 

+ \ a(l ll l m T lrn + 7 "« 2 T 00 ) 
A* = ^^ kl r p ^(T ml T pkj - T nij T pkl ) - 2^7r'%« - ^I7 2 « 2 r 00 , (54) 
which can be written in a simpler way as 

Aa = -±9, ia ,i - , + ^(T u + *a*T™) + \hj?a, 3 - \f^9 M , (55) 



A* = . _ 1 ^ ^VT 00 . (56) 

4* ' ' 2* 3 2 



Here, we have used the fact that the noncovariant trace of the Ricci tensor at this order is given just by, 

R* = ^^ + ^*« + ^*,«- (57) 

Finally, the set of CFC+ equations is closed by the elliptic equation for the transverse and traceless part of the 
spatial metric, 

AhJT = Bij - - V\i - V* d + ^S tJ V m , m , (58) 

Al ' • = />',,, [ - !>>.,. (59) 

where B = Bn and now the following definition for B^ should be applied, 

= f.jttj + >!'..,(»., - uTij. (60) 
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In order to arrive to this expression, we have integrated by parts terms with second derivatives of the lapse and the 
conformal factor, taking into account that the full-derivative terms are eliminated when computing the transverse 
and tracefree part, as it is done in the right-hand side of Eq. (|58]) . These CFC+ equations have been compared to 
those presented in Ref. [l2[ and obtained exact agreement. For such comparison, one has to consider that, as already 
shown in the previous section, = 1 + 2U/c 2 and a = 1 — U/c 2 for the 2PN terms beyond CFC. 

In conclusion, we have shown that our proposed 4PN-exact formalism covers naturally both CFC and CFC+ 
approximations. Regarding the 2PN-exact approach presented in Sec. IIV1 it also coincides with the CFC and CFC+ 
cases up to the mentioned PN-level. The key difference between the CFC+ scheme and the one of Sec. IIVI is that 
whereas we have kept all the nonlinear terms present in all the equations, in CFC+ approach one just keeps those 
terms that are of the corresponding (2PN) order. Even though, our 2PN elliptic framework of Sec. IIVI can not "fully" 
cover the CFC case since we have assumed that 7r^ T = whereas in the standard CFC scheme the definition (f50|) is 
used. Only under spherical symmetry conditions both hJ T and tt^ t vanish. 



VII. A SPECIFIC MATTER MODEL 

In this section we specify the matter model to the particular case of a barotropic perfect fluid [in this section, c = 1]. 
Therefore, the stress-energy tensor will be given in terms of the fluid four- velocity u 11 , pressure p, proper mass density 
p, and specific enthalpy h in the following way, 

T^ v = p(l + K)u^u v + pg^, (61) 

where the conservation law V M (pw M ) = holds with V M the four-dimensional covariant derivative. 

We need the following components of the stress-energy tensor density {—g) x l 2 T^ v that appear in the equations, 

a7 i/2 T oo = p n +h)u t_f^P j (62) 

a 

aj i/2 T o. = p ^i + h) Uj , (63) 

ai x l 2 T l3 = p^l + h)u i l ^- t +a 1 1 l 2 p ll3 , (64) 

where the four-velocity is decomposed as u M = (u*, u 1 ), = (u t , Ui) and the definition p» = a^^v^p has been made. 
From the normalization of the four- velocity u^u^ = —1, it is easy to see that it* is given by u* = (1 + j^u l u 3 )^ 2 /a. 

In the case of a barotropic perfect fluid, where dp — p dh holds, the equation of state p = p{p) and the conservation 
of stress-energy tensor, 

V M T^ = 0, (65) 

give rise to all the equations of motion for the matter. The independent equations of motions can be cast into the 
form [Hill, 

dtp* = -di(p*v l ), (66) 

d t P t = -djiPitf) - a 4 (a 7 1/2 p) + ^l 1/2 T^d l9 ^, (67) 

where v l = u l ju l and Pi = p*(l + h)ui. Introducing Wi = Pi/p*, i.e., Wi = (1 + h)ui, the latter equation of motion 
can be written 

1 1/9 aj 1/2 Ti- ll/ 
d t w l = -v 3 d j w l d t (a'f L/z p) H d.g^. (68) 



With the aid of the relation 

al " W ' rmS. « 



[(1 + h) 2 +j i iw l w J } 1 / 2 
all matter variables can be reduced to the independent ones (p*,Wi) or (p*,Pi). 
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VIII. CONCLUSIONS 

In this paper we have developed a 4PN-exact approximation to the field equations of General Relativity which turns 
out to be fully exact in the conformal flat case and 5PN exact for stationary configurations. The elliptic equations for 
the lapse a, the shift f3 l , and the conformal factor \fr are exact (just the well-known equations from the ADM formalism 
but in some more explicit form): Eqs. (fl~7|) . ([TO . and (|2Tj) . respectively. Only the wave equation for hj^ (|47f and 

the definition for its conjugate momentum tt^ t (|42]l are approximated. In order to obtain the mentioned hyperbolic 
equation, an implicit transition from a Hamiltonian to a Routhian framework, regarding the independent gravitational 
degrees of freedom, has taken place Q . This equation has to be solved under the condition of no incoming radiation. 
The construction of various transverse-traceless objects resulted in several auxiliary vectors tt 1 , V z , and W l obeying 
the Eqs. (jUJ), (|33j), and flU), respectively. 

Our proposal is to solve the commented set of equations iteratively. In particular, the first step is to solve the 
elliptic equations and the equations of motion of the matter for all times prior to the initial value slice imposing 
h% T = and irj^ = 0. Then, using this result, the hyperbolic equation for h^ T as well as the elliptic equations again 
are being solved which gives rise to a 2[2.5]PN-exact approximation to the Einstein field equations. The obtained 
2PN-exact solution can be used to correctly compute 7r^ T through its truncated definition (|42l) up to 3PN order. 
With this information at hand, the sources of the mentioned system of equations can be calculated to one PN-order 
of precision higher. This fact permits, by solving again the system, to obtain a 3[3.5]PN-exact solution. This process 
can be repeated up to obtain a 4[4.5]PN-exact solution. The sources have been manipulated so that they decay fast 
enough (~ 1/r 4 ) at asymptotic infinity. Therefore the system of equations we present is well suited for an eventual 
numerical evolution. At any instant of time, the obtained solutions can be used as initial data to be evolved according 
to the full Einstein equations or, the solutions can be treated as approximate solutions of the full Einstein equations 
throughout all times. 

The presented approach covers recently developed and explored approaches including outgoing radiation [l3l - [T5l . |27| . 
For binary compact objects, metric coefficients for the coordinate conditions of the present paper are known in closed- 
analytic form through 2[2.5]PN order in the near zone [1, H, H, HI] . Those have already been used for further iterations 
in, e.g., 0, [H. 

Finally, we have also analyzed the 2PN-exact waveless approximation to General Relativity that is obtained by 
imposing the assumptions = 0. These conditions are quite interesting since they can be considered as 

the conjugate assumptions to the well-known conformal flat conditions (hj^ = hJ } T = 0). A maximum elliptic system 
without any partial time derivatives of the gravitational field variables that one can obtain within General Relativity 
is dropping hj^ in Eq. (|3U)l and 7r^ Tj0 in Eq. (|3"0j) . Hereof, the simplified waveless approach suggested in Sec. [V] 

results by just droppin g hj ? in Eq. f|42[) . A further truncation is obtained by putting ir^ T = which is the waveless 
approximation of Sec. IIVI All these waveless approximations are 2PN exact where the latter one is the simplest 
of them. In order to obtain more accurate truncated versions of General Relativity, hyperbolic equations should be 
considered, as it is done in the above commented 4PN-exact approach. 
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